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Quantum Monte Carlo (QMC) and density-matrix renormalization group (DMRG) methods 
are used to study the coupled spin-pseudospin Hamiltonian in one-dimension (ID) that models 
the charge-ordering instability of the anisotropic Hubbard ladder at quarter filling. We calculate 
the temperature dependence of the uniform spin susceptibility and specific heat as well as the 
spin and charge excitation spectra of the system. We show that there is a parameter and 
temperature region where the spin degrees of freedom are separated from the charge degrees 
of freedom and behave like a ID antiferromagnetic Heisenberg model, and that, outside this 
parameter region and above a crossover temperature, the spin excitations are largely affected 
by the charge fluctuation. We argue that observed anomalous spin dynamics in the disorder 
phase of a typical charge-ordered material a'-NaV2 05 may possibly be a consequence of this 
type of spin-charge coupling. 
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1. Introduction 

Charge-ordering (CO) instability has recently been one 
of the major topics in the field of strongly correlated elec- 
tron systems. Here, elucidation of the observed anoma- 
lous behaviors of electrons associated with the CO phase 
transition has been the central issue. This includes ques- 
tions on the charge dynamics above the transition tem- 
perature Tco as well as on the CO spatial patterns real- 
ized below Tco- A well-known example is the vanadate 
bronze Q;'-NaV205 where the system may be modeled as 
a lattice of coupled ladders (or a trellis lattice) at quarter 
filling. Strong intersite Coulomb interaction between 
electrons is believed to be the origin of the CO instabil- 
ity.^' In this material, the CO with a zigzag ordering 
pattern is observed below Tco — 34 K,^"^"^ and associ- 
ated with this, a number of anomalous behaviors, which 
can be related to the slow dynamics of charge carriers (or 
charge fluctuation), have been observed above Tco.^°"^^^ 
Anomalous response of the spin degrees of freedom has 
also been noticed.^' i^^^i) seems therefore quite natural 
to wonder how in such systems the spin degrees of free- 
dom behave near the CO phase transition when they are 
on the slowly fluctuating charge carriers. In this paper, 
we thus consider the issue: what are the consequences 
of charge fluctuation at T > Tco to the spin degrees of 
freedom? 

One of the simplest models that allow for such situa- 
tion is the anisotropic Hubbard ladders at quarter filling 
with the strong intersite Coulomb repulsion. Wc here use 
an effective Hamiltonian written in terms of the spin and 
pseudospin (representing charge degrees of freedom) op- 
erators.^' ^^'^^'^'^-' This Hamiltonian is derived from the 
Hubbard ladder model by the perturbation theory^' 
where the hopping parameter between the rungs of the 
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ladder is assumed to be small compared with the onsite 
and intersite Coulomb repulsions as well as the hopping 
parameter in the rung (i.e., the anisotropic ladder). Al- 
though the CO is not realized in this model (since it is 
the ID quantum-spin model), we can simulate anoma- 
lous behaviors of the spin degrees of freedom under the 
strong charge fluctuation. We will apply the quantum 
Monte Carlo (QMC) method to this model to calculate 
the temperature dependence of the uniform spin sus- 
ceptibility and the spin and charge excitation spectra, 
thereby clarifying consequences of the interplay between 
its spin and charge degrees of freedom. The density- 
matrix renormalization group (DMRG) method with the 
finite-temperature algorithm will also be used to calcu- 
late the temperature dependence of the specific heat of 
the model. 

We note that, in this coupled spin-pseudospin model, 
the spin exchange interaction is necessarily associated 
with the charge excitation; i.e., the spin excitations can- 
not occur without making the exchange of the pseu- 
dospins. We will then show that nevertheless there is a 
parameter and temperature region where the spin de- 
grees of freedom behave like a ID antiferromagnetic 
Heisenberg model; i.e., the spin degrees of freedom are 
'separated' from the charge degrees of freedom in this 
region. We will moreover show that the spin system be- 
haves in different manner depending on whether the tem- 
perature T is below or above a crossover temperature T* 
that is related to the pseudospin excitations; at T < T*, 
it behaves like a ID antiferromagnetic Heisenberg model 
with a T-independent effective exchange coupling con- 
stant Jeff with large renormalization, whereas at T > T*, 
Jeff decreases rapidly with increasing T, where the ef- 
fective Heisenberg-model description ceases to be valid. 
Because the parameter values for a'-NaV205 are outside 
the region where the spin-charge separation is complete, 
we will argue that some experimental data may possibly 
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be interpreted as consequences of this type of spin-charge 
coupling. 

This paper is organized as follows. In §2, we define the 
coupled spin-pseudospin model that describes the spin 
and charge degrees of freedom of the anisotropic Hubbard 
ladder at quarter filling. Some details of the method of 
calculation are also given. In §3, we present the results 
of calculation which include the staggered susceptibil- 
ity for pseudospins, the spin and pseudospin excitation 
spectra, and the temperature dependence of the uniform 
spin susceptibility and specific heat. Discussion on the 
experimental relevance to a'-NaV205 and summary of 
the paper will be given in §4. 

2. Model and Method 

Our effective spin-pseudospin Hamiltonian may be 
written as a sum 



H = TIq + HsT 



(1) 



of the quantum Ising Hamiltonian for pseudospins 

L L 

Wo = Ji ( - f E + E ^i^i+i) (2) 

i=l i=l 

and the spin-pseudospin coupling term 

WsT = J2 E i^i ■ Si+i - 4) {T^Tr^x + H.c.) . (3) 

1=1 

The standard notation is used here. Sj and Tj are, re- 
spectively, the spin and pseudospin operators of spin- 1/2 
at site i, where = —1/2 (+1/2) means the electron is 
on the left (right) site on the rung of the ladder. L is the 
system size and periodic boundary condition is assumed. 
Ji is the energy scale of the pseudospin system and J2 is 
the coupling strength between the spin and pseudospin 
systems. 

From the second-order perturbation theory,^' ^^■^■^^ we 
have the relations Ji = 2V\\ and J2 = ^^/V±, where 
t|| and V|| {t± and V±) are the nearest-neighbor hopping 
parameter and Coulomb repulsion of the leg (rung) of the 
ladder, respectively. We should then have Ji > J2, which 
we assume throughout the present work. We also assume 
the onsite Coulomb repulsion to U 00. Relative 
strength of the transverse field applied to the pseudospins 
is measured hy g = 4i_\_/J\ = 2t±/V\\. Note that g in the 
quantum Ising model represents the relative strength of 
the fiuctuation of a charge in the rung: if we assume 
one electron in a rung, we have the prefactor gJ\/2 in 
the first term of eq. (2), which is the difference between 
the energies of the bonding and antibonding levels of the 
rung, 2t_i_. Thus, if g (or t±_) is large the electron is stable 
in the bonding level of the nmg, but if g (or t^_) is small 
the effect of Vn easily leads the system to CO. 

We use the conventional world-line QMC method for 
the analysis of the model. We use a 32-sitc cluster (where 
a site contains a spin and a pseudospin) with periodic 
boundary condition; the cluster-size dependence of the 
calculated results are examined by using clusters up to 
96 sites but we find no significant size dependence in 
the results. Because the model does not conserve the to- 
tal pseudospin, we have examined a number of ways of 
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Fig. 1. Tbiupcraturc dependence of the staggered susceptibility 
for pseudospins XT ("■) calculated for the coupled spin-pseudospin 
Hamiltonian. 



the spin fiips and confirmed that available analytical re- 
sults are reproduced correctly. The maximum-entropy 
method is used to calctilate the dynamical quantities like 
the spin and pseudospin excitation spectra. The DMRG 
method with finite-temperature algorithm^^' is also 
used for the calculation of the temperature dependence 
of the specific heat of our model; the method enables us 
to access to lower temperatures. 

3. Calculated Results 

3. 1 Staggered susceptibility for pseudospins 
The response function is defined as 



Xi] 



d\{{Sn-i\)St)-{S^,){St)) (4) 



where Sj{—i\) is the Heisenberg representation of 
and (• • • ) is the canonical average. Xij is Fourier trans- 
formed to the g-dependent susceptibility x('z)i which we 
calculate by the QMC method; the g — » limit gives the 
uniform spin susceptibility x(^) and the staggered sus- 
ceptibility is defined as x(g) at g = tt. In the following, we 
use the subscripts S and T as in xs('7) and Xt(9), which 
stand for the spin and pseudospin degrees of freedom, 
respectively. 

The phase diagram of the quantum Ising model Hq 
is well known;^^) at T = there is a long-range order 
for 5 < 1 (51 = 1 is a quantum critical point), which 
corresponds to the zigzag (or 'antiferromagnetic') CO. 
The calculated staggered susceptibility for pseudospins 
is shown in Fig. 1, where we find that it shows divergent 
behavior at T ^ for g < 1. The dispersion relation 
of the pseudospin excitation observed in the calculated 
dynamical structure factor (shown in Fig. 2, see below) 
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agrees well with the exact result 

Ji 



9=2.0 



g=2.0 



= yVl + ff^ + 2.9COS(7. 



(5) 



We find in Fig. 1 that the inclusion of the coupling term 
HsT, which introduces the quantum fluctuation via the 
factor T^T~ , siipprcsscs the divergence. Thus, wc may 
say that the inclusion of the spin degrees of freedom in 
the quantum Ising model for pseudospins leads to the 
unstable long-range order of pseudospins. 

3.2 Spin and pseudospin excitation spectra 

The dynamical pseudospin structure factor 5t(5, w) is 
defined as 

ST{q,r) = ^Ee-*«('-^-'--)(T;^(r)T;^(0)) (6) 

ij 

1 

ST{q,T) = - dui ST{q,co)K{u!,T) (7) 

Jo 

K{lj, t) = e-'^'^ + e-'^^'^-^) (8) 

where S'T(g, t) is the Fourier transform of the imaginary- 
time correlation function. We use the maximum entropy 
method for the inverse Laplace transformation (or an- 
alytical continuation) to obtain S'liq,'-^) from S't(9,t). 
The dynamical spin structure factor Ss{q,io) is similarly 
defined by replacing the pseudospin operator with the 
spin operator S^. 

The calculated results for the pseudospin excitation 
spectra at low temperature {ksT = O.IJ2) are shown in 
Fig. 2, where wc find that the spectra arc under strong 
infiuence of the spin-pseudospin coupling term J2- With 
increasing the coupling strength J2/J1, the peak of the 
pseudospin spectra shifts to higher energies and simul- 
taneously the spectra are broadened. Thus, the lower- 
energy edge of the peak is not affected strongly by the 
coupling strength J2, at least when g is large. It seems 
reasonable to suppose that the scattering of the pseu- 
dospin excitations due to spin excitations causes the 
broadening of the spectra. 

The calculated results for the spin excitation spectra at 
low temperature are shown in Fig. 3, where we find that, 
in contrast to the pseudospin spectra, the spin excitation 
spectra change very little; i.e., the peak position, width, 
as well as the shape of the spectra are not affected by the 
parameter Ji when .9^1. When g is small, however, the 
peak position is slightly shifted to lower energies with 
increasing the value of Ji (see Fig. 3 (a)). 

The dispersion relation of the spin and pseudospin ex- 
citations calculated at low temperature are summarized 
in Fig. 4, which are obtained as the momentum depen- 
dence of the peak position of the spectra. For compari- 
son, we show the dispersion of the quantum Ising model 
in Figs. 4 (a) and (c); the gap opens when g > 1, which 
is closed at q = n when g ^ 1, leading to the 'antifer- 
romagnetic' long-range order (or zigzag CO).^^'^^-* We 
note that the gap remains open irrespective of the value 
of g when wc include the coiipling term J2. In the left 
panels and right panels of Fig. 4, we present the same 
dispersion relations Wq, but in a different energy scales, 
i.e., u)q/Ji and u)q/J2. We find that, unless g is small, the 
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Fig. 4. Dispersion relations of the spin (open symbols) and pseu- 
dospin (solid symbols) excitations calculated at ksT = O.IJ2. 
Note that the same data at g = 2 (at g = 1) are plotted in 
(a) and (b) (in (c) and (d)) in different energy scales Ji and 
■J2. The dotted line in (a) and (c) is the dispersion relation for 
the quantum Ising model cq. (5), and that in (b) and (d) is the 
scaled dispersion relation for the ID antiferromagnetic Heisen- 
berg model eq. (9). 



spin excitation spectra are always inside the charge gap, 
i.e., inside the gap of the pseudospin excitation spectrum. 
Thus, when the charge gap is large, the energy scale of 
the spin excitations is separated from the high-energy 
charge excitations. With decreasing g, however, the en- 
ergy of the charge excitation decreases at the momentum 
q = TV to couple with the spin excitations. 

In Fig. 4 (b) and (d), we find that, for .g > 1, the 
dispersion of the spin excitation spectra scales very well 
with J2; i.e., it does not depend on the value of J\. More 
quantitatively, the dispersion of the calculated spin ex- 
citation spectra is fitted well with the dispersion of the 
ID antiferromagnetic Heisenberg model 



Wg/ J2 = 0.6 X — sin g 



(9) 



if we include the factor 0.6 as in eq. (9). The factor is 
independent of Ji for > 1 and at low T. 

These results suggest that at low temperatures there 
is a parameter region where the spin degrees of freedom 
behaves independently from the pseudospin degrees of 
freedom; it is when 5 ^ 1 and the gap of the pseudospin 
excitation spectra is large, inside of which there is a spin 
excitation spectra. Thus, we suggest the validity of the 
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Fig. 2. Dynamical pseudospin structure factor ST{q,oj) for the coupled spin-pseudospin model calculated at ksT = 0.1^2- The results 
at J2/J1 = are for the quantum Ising model. The peaJj at u) = for J2/J1 > in the uppermost panel of (a) and (b) is spurious, 
which is due to the error of the maximum entropy method. 
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Fig. 3. Dynamical spin structure factor S'g(g,t<j) for the coupled spin-pseudospin model calculated at ksT = O.IJ2. 



decoupling of the coupling term of the Hamiltonian as 
HsT ^ J2 E (r^^^i + H-c.) (S, • S,+i - -) (10) 

i=l 

with 

{T,+T-^,+ii-c.)c^0.6 (11) 

which leads to the effective Heisenberg-model descrip- 
tion of the spin degrees of freedom of our model. Here, it 
should be noted that in gcncr al the factor < T+T-^_^ + 
H.c. > at zero temperature takes the value 1 when g ^ 
and 1/2 when g ^ oo; thus the above value 0.6 reflects 
the effect of quantum fluctuations of pseudospins, which 
is strong already at g > 1. It should also be noted that 
with increasing temperature the value of the factor de- 
creases to (due to thermal fluctuations), either very 
rapidly when g is small or rather slowly when g is large. 



We flnd that the essential features of < T^T^^^ +H.c. > 
are contained already in a two-site (or dimer) model of 
eq.(l), a minimum model reflecting the spin-pseudospin 
coupling. This is evidently shown in Fig. 5. 

3.3 Uniform spin susceptibility 

To see the validity of the effective Heisenberg-model 
description further, in particular for its temperature de- 
pendence, we calculate the temperature dependence of 
the uniform spin susceptibility for the coupled spin- 
pseudospin Hamiltonian. The results are shown in Fig. 6, 
where comparisons are made with the uniform suscep- 
tibility for the system of free spins and with that for 
the ID antiferromagnetic Heisenberg model. We find 
that the temperature k-QT/J^ at which J2Xs{T) shows 
a maximum is lower than that of the ID antiferromag- 
netic Heisenberg model; it becomes lower with decreasing 
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Fig. 6. Temperature dependence of the uniform spin susceptibil- 
ity xs calculated for the coupled spin-pseudospin Hamiltonian. 
The solid and dotted curves are the uniform susceptibility for 
the system of noninteracting 5 = 1/2 spins and that for the ID 
antiferromagnetic Heisenberg model, respectively. 



the value of g or with increasing the value of J1/J2. In 
other words, the deviation from the Heisenberg model is 
large when the quantum fluctuation of the pseudospins 
is small, which occurs when g is small or Ji is large. 

Now, let us analyze the data more precisely. In or- 
der to do this, we fit the results with the temperature 
dependence of the spin susceptibility of the ID antiferro- 
magnetic Heisenberg model, the so-called Bonner-Fisher 
curve;^^-' i.e., we introduce the T-dependent effective ex- 
change coupling constant Joff(T) and we determine the 
values so as to fit the calculated uniform spin suscepti- 
bility xs(r). If the values of Jes thus obtained do not 
depend on T, it follows that the spin degrees freedom of 
our spin-pseudospin model is reduced to a ID Heisenberg 
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Fig. 7. Effective exchange coupling constant Jeff(T) estimated 
from the fitting of the calculated uniform spin susceptibility to 
the Bonner-Fisher curve. Note that the same data are plotted 
as a function of k^T/Ji (left panels) and of k^T/ J2 (right pan- 
els), whereby a scaling behavior is seen in the latter. The arrows 
indicate the crossover temperature T*. The solid lines are the 
guide to the eye. 



model 



Wspin — "^off ^ ^ Si ■ S 



(12) 



at least for the response to the uniform magnetic field. 
The results are shown in Fig. 7. We find that the esti- 
mated value of Joff(r) is indeed a constant for temper- 
atures below k-QT < 0.7 J2 at g = 2. A crossover tem- 
perature T* (= O.7J2) is thereby defined. The effective 
exchange coupling constant thus deduced takes a value 



Joff ~ O.6J2 



(13) 



at r < T*; this value is consistent with the value es- 
timated from the dispersion relation of the spin exci- 
tation spectra (see §3.2). We find that also at g = 4 
the scaling behavior holds up to a higher temperature 
{^bT < O.8J2), but with a slightly smaller value of Joft 
(see Fig. 7 (d)), demonstrating the vahdity of the effec- 
tive Heisenberg- model description at T < T*. At g = 1, 
however, the temperature region where Jeff (2^) takes a 
constant value is already very small, although the value 
is still Jeff ~ 0.6 J2 at T ~ K, and at g = 0.5, the value 
of Jeff at T ~ K deviates largely from Jes = O.6J2 (or 
decreases strongly when J1/J2 is large), where the effec- 
tive Heisenberg- model description completely fails. We 
thus find that T* thus deduced is insensitive to Ji , scales 
well with J2, and depends strongly on g (i.e., T* 00 
at (7 ^ 00 and T* ^ at g ^ 1). Note that if we assume 
eq. (10) the behavior should come from the pseudospin 
fluctuation {T^Tfl^^ 4- H.c), and actually we find that 
very rough tendency in the parameter and temperature 
dependence is seen in the results for the dimer model (see 
Fig. 5) although the scaling behavior and the presence 
of T* are not seen. 
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We note here that the crossover temperature T* 
roughly scales with J2 rather than Ji, as seen in Fig. 7. 
One might suppose that it should scale with the size of 
the charge gap: i.e., up to temperatures corresponding to 
the energy of the lowest charge excitations, with which 
the pseudospins can excite, the spin excitations may be 
written in terms of the ID antiferromagnetic Heisenberg 
model. However, as we have discussed in §3.2, the size of 
the charge gap (if it is defined as a low-energy edge of 
the peak) shows a rather complicated behavior and does 
not simply scale with either J2 or Ji. The naive picture 
thus does not hold. It may be said however that, since 
there is no other excitations available, the deviation from 
the ID Heisenberg-model description is necessarily due 
to the pseudospin excitations. 

3.4 Specific heat 

Finally, we present the calculated results for the tem- 
perature dependence of the specific heat C by the finite- 
temperature DMRG method. ^^'^^^ The results are shown 
in Fig. 8, where we find that the curves have a two-peak 
structure: e.g., a,t g — 2 and J1/J2 = 4, a rather sharp 
peak appears at a low-temperature region k^T ~ O.3J2 
(which scales with J2) and a broad peak structure ap- 
pears at a high-temperature region ksT > 0.8 J2 (which 
scales with Ji), each of which corresponds to the spin 
and pseudospin excitations, respectively. We also find 
that the shape of the low-temperature peak can be fitted 
very well with the temperature dependence of the specific 
heat of the ID antiferromagnetic Heisenberg model with 
the effective exchange coupling constant J^g ~ O.6J2 at 
(7 = 2; this value is in accord with the value estimated 
in §3.2. The temperature at which the deviation in the 
fitting occurs is at k-^T / J2 — 0.5 — 1 depending on the 
value of (7, which is also consistent with the estmate from 
the temperature dependence of the uniform spin suscep- 
tibility. The results thus demonstrate the separation of 
the energy scales between spin and pseudospin degrees 
of freedom and the presence of low-energy 'magnetic' en- 
ergy scale as has been pointed out in ref.^^-* 

4. Summary and Discussion 

We have calculated the spin and pseudospin excitation 
spectra and the temperature dependence of the uniform 
spin susceptibility of the coupled spin-pseudospin Hamil- 
tonian by using the QMC method. We have also calcu- 
lated the temperature dependence of the specific heat 
of the model by the finite-temperature DMRG method. 
We have first shown that, when the pseudospin quan- 
tum fiuctuation is large {g > 1), the dispersion relation 
of the spin exitation spectra of our model at low tem- 
peratures agrees well with that of the ID antiferromag- 
netic Heisenberg model with the renormalized effective 
exchange coupling constant Jeff — O.6J2 that is inde- 
pendent of the energy scale of the pseudospin system 
Ji. Here, the spin excitation spectra is well inside the 
charge gap, and thus the spin degrees of freedom are sep- 
arated from the charge degrees of freedom. We then have 
shown that the temperature dependence of the uniform 
spin susceptibility of our model is well described again 
by the ID antiferromagnetic Heisenberg model with the 




Fig. 8. Temperature dependence of the specific heat C calculated 
for our coupled spin-pseudospin Hamiltonian. 



same effective exchange coupling constant Jeff — O.6J2. 
This description is valid up to the crossover tempera- 
ture T* that is related to the pseudospin excitations of 
the system and roughly scales with J2 unless the quan- 
tum fluctuation of the pseudospins is small (5^1). We 
have shown the appearence of the two-peak structure 
in the temperature dependence of the specific heat and 
have confirmed the presence of the low-energy magnetic 
energy scale. We have thus demonstrated the validity 
of the effective Heisenberg-model description of the cou- 
pled spin-pseudospin model for the quarter-filled ladders. 
It then follows that the coupling between the spin and 
pseudospin degrees of freedom, which occurs at g < 1, 
leads to the anomalous spin and charge dynamics of the 
system. 

Finally, let us discuss some possible experimental rel- 
evance of our results. We have rough estimates of the 
values of the physical parameters appropriate for a'- 
NaV205; according to ref.,'^^ we have t\\ ^ 0.14 eV, 
ti_ ~ 0.30 eV, and I^, - VI ~ 0.8 eV, which lead to 
Ji - 1.6 eV, J2 - 0.10 eV, and g ~ 0.75. We thus 
find that the real material may be in the region of g < 1, 
where the spin degrees of freedom are not completely sep- 
arated from the charge degrees of freedom. The anoma- 
lous response of the spin degrees of freedom may there- 
fore be expected. We here point out, e.g., that the value 
of Jeff estimated from the uniform susceptibility observed 
in experiment (which takes the value ^ 600 — 700 K 
at T K) indeed decreases with increasing tempera- 
ture,^^ which is consistent with the results of our calcula- 
tion. Our explanation is that the effective spin exchange 
coupling constant depends primarily on the fiuctuation 
of pseudospins, the temperature dependence of which is 
strong above the crossover temperature (T* ~ K in this 
material) where the spin-charge coupling becomes rele- 
vant. The reported^*^-* temperature dependence of the nu- 
clear spin-lattice relaxation rate 1/Ti in this material is 
also interesting in the present respect. The measured^^^ 
strong temperature dependence of the integrated optical 
spectral weight have also been discussed in terms of the 
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destruction of short-range spin correlations which oc;curs 
as the temperature is increased. Although the real ma- 
terial a'-NaV205 is modeled well as a 2D trellis-lattice 
system in particular for its charge degrees of freedom, we 
here want to point out that observed anomalous spin dy- 
namics in the disorder phase of this typical CO material 
may possibly be a consequence of this type of spin-charge 
coupling. 

Because the anomalous charge dynamics has been no- 
ticed also in other transition-metal oxides^^^ and some 
organic systems, we hope that the present study will 
stimulate further researches on the intriguing interplay 
between the spin and charge degrees of freedom of 
strongly correlated electron systems with CO instability. 
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